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C$ ■ Abstract 

^p 1 i We study the equilibrium phase diagram of a generalized ABC model on 

■ an interval of the one-dimensional lattice: each site i = 1, . . . , N is occupied 
Ch I by a particle of type a = A, B, C, with the average density of each particle 

species N a /N = r a fixed. These particles interact via a mean field non- 

1 1 1 reflection-symmetric pair interaction. The interaction need not be invariant 

under cyclic permutation of the particle species as in the standard ABC 

^2 ! model studied earlier. We prove in some cases and conjecture in others that 

£^ ■ the scaled infinite system N —> oo, i/N — > x G [0, 1] has a unique density 

| profile p a {x) except for some special values of the r a for which the system 

■ undergoes a second order phase transition from a uniform to a nonuniform 
| periodic profile at a critical temperature T c = "Sy/rAtBTc /27T. 

\6 \ 
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1 Introduction 

X 

d ' The standard ABC model on an interval was considered in [I]. It is an equilibrium 

system on a ID lattice of N sites with closed boundary conditions. Each site is 
occupied by one of three types of particles, denoted by A, B, and C. The particles 
interact via a cyclic mean field type pair potential which is however not spatially 
reflection symmetric. In this paper we generalize this model by introducing an 
additional interaction in which each of the particle types moves in a separate 
background potential that depends linearly on position. This breaks the cyclic 
symmetry in the standard ABC model. The equilibrium state of the standard 
ABC model may also be obtained as the steady state for certain nearest neighbor 
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exchange dynamics [U US E] , and the generalized model considered here may be 
obtained by a modification of the exchange rates; see Appendix [X] for details. 

To define the model we introduce the occupation variables rj a (i), with r) a (i) = 
1 (0) if site i is (is not) occupied by a particle of type a. As each site is occupied 
by exactly one particle, 

!>«(*) = i. (i.i) 



The energy of a configuration ij is defined to be 

1 N I N \ 

a i=l \j=l J 

Here a + 1 corresponds to the species following a in the ABC cyclic order, 0(j — 
i) = 1 (0) for j > % (j < i), and the £ a may be thought of as constant background 
electric fields, with £ a acting on particles of type a. We will consider the model 
in the canonical ensemble, with specified particle numbers of each type: 

N 

$>«(») = #«, ]TiV Q = iV. (1.3) 

i=l a 

The canonical Gibbs measure for this system is then given by 

MmN a ) = ^e-^), (1.4) 

with Z the usual canonical partition function. We will assume N a > for all a 
throughout; if one of the particle species is absent, the model simply reduces to 
the weakly asymmetric simple exclusion process (WASEP) [HO2]. 

Note that if one adds the same constant to each of the £ Q then the energy 
(jl.2p is only changed by an overall constant. We may therefore set ^2 a Ca = 0, 
without loss of generality. We will refer to the case where £4 = £b = £c = as 
the "standard" ABC model. This is the model considered in pQ. 

The energy (jl.2p may also be written in a different form, in which the contribu- 
tion of the external fields is expressed through a modified mean field interaction. 
Using (jl.ip and fjl .3j) we have 

N N N 

i=l i=l j=l (1.5) 

Substituting (jl.5p in (jl.2p and rearranging sums, we obtain 

N N 

N 



1 ( N N 



i=1 i =1 (1.6) 

+ ta^(N + N a+2 -N a+1 + l)Y 
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Figure 1: Different regions of the parameter space of the v a , plotted in the va-vb 
plane. 

where the v a are given by 

Va = ^ U + - £a+2\, (1-7) 

so that 

£>« = !■ (1-8) 

a 

The term ^ a £ Q N a (N + iV Q+ 2 — N a+ i + 1) /2 is independent of the configura- 
tion, and may be ignored in the canonical ensemble with the N a fixed. 

It will be convenient to consider the fundamental parameters of the model, 
in addition to the particle numbers N a of each species, to be the v a rather than 
the fields £ Q , as the v a are more directly related to the physical behavior of the 
model. We divide the space of these parameters into three regions; see Figure HJ 
plotted in terms of va and vb, as these fix vc by (jl.8|) . In region I, v a > for all 
a. In regions II and III, v a < for one or two values of a, respectively. In the 
standard model va = vb = vc = 1/3. 

We see from (II. 6ft that each of the v a determines the energetically preferred 
alignment of the two other particle species adrl. Effectively there is a contribution 
of 3 v a /N to the energy every time any pair of particles of species a + 1 and a + 2 
are not cyclically aligned, that is, whenever a particle of type a + 2 precedes one 
of type a + 1. If for example va > 0, a pair of B and C particles will have a lower 
energy arranged as B ■ ■ ■ C than as C ■ ■ ■ B. If va < the preferred arrangement 
is reversed, and the configuration C ■ ■ ■ B will have a lower energy than the usual 
cyclic ordering B ■ ■ ■ C. This determines the ground states of the system, when 
/3 — > oo, as described in Appendix [Al 
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Background 

The standard ABC model was introduced by Evans et al [21 [3] and in the 
form we use by Clincy et al [BJ. This model was originally considered on the ring, 
i.e. with the boundary conditions periodic rather than closed, by specifying a 
dynamics consisting of asymmetric nearest neighbor exchanges between particles 
of different species. The stationary state of this dynamical system on the ring 
is generally not an equilibrium one. Its properties have been studied extensively 
in [3] and in [7J [HJ [TOJ CD] ■ In the special case that the number of particles 
of each species is the same the stationary state of the model defined on the ring 
is a canonical Gibbs measure with the energy given by (jl.6p . with v a = 1/3 
for all a. The stationary state of the dynamical model defined on the interval, 
or equivalently on a ring with exchanges across one bond blocked, is always the 
canonical Gibbs measure regardless of the number of particles of each species and 
for equal N a is identical to that on the ring. 

The equilibrium properties of the standard ABC model on the interval were 
obtained exactly in [1] , see also [12] . In particular it was shown there that (in the 
thermodynamic limit) the system has a unique state (density profiles) whenever 
the average densities N a /N = r a are not all equal to 1/3, i.e. r a ^ v a for some a. 
When r a = 1/3 the system undergoes a second order transition at j3 = j3 c = 2tt\/3 
from a uniform density profile to a periodic profile [8] . There are thus for (3 > (3 C 
a continuum of phases (density profiles) specified by a rotation angle 0. The 
results we derive here for the case when the v a are not all equal to 1/3 are 
more restricted. They suggest however that the phase diagram for general v a is 
qualitatively similar to that of the standard ABC model, with a phase transition 
only at r a = v a . 

The outline of the rest of the paper is as follows. In Section [2] we discuss 
general properties of the macroscopic system in the scaling limit. Section [3] de- 
scribes basic properties of the solutions of the Euler-Lagrange equations. Special 
cases are considered in Sections U] and for different values of the parameters. 
In Appendix [X] we discuss properties of the microscopic model specified by f) 1 . 2 j) . 
Appendix [B] discusses the connection between the generalized ABC model and 
the Lotka-Volterra family of ODE systems, and in Appendix Owe discuss some 
restrictions on solutions of the Euler-Lagrange equations for the case r a ^ v a for 
some a, with the v a in region I. 

2 Scaling limit 

The main goal of this paper is to study the phase diagram of the equilibrium 
system with energy (jl.2p when N becomes macroscopic. For this we consider the 
scaling limit in which 

iV^oo, N a /N^r a , i/N -»• x £ [0, 1], (2.1) 

so that r a is the average density of the particle species a. In this limit the state 
of the system is described by density profiles p a (x), a = A,B,C, where p a (x) 
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represents the density of particles of type a at a position x. These density profiles 
satisfy the constraints 

< p a (x) < 1, / dx p a (x) = r a , y2p a (x) = l. (2.2) 

Let 5({p}) and £({p}) be the entropy per site and energy per site in the 
scaling limit. As in [T], see also [SI EE], these are given by 

$({p}) = Y] [ dx p a (x)logp a (x), (2.3) 
£{{p\) = y2 dx d v e (y ~ x )Pa{x)Pa+2{y) 

a JO JO 

+ y / dx x£ a p a (x) + constant, (2-4) 

a JO 

= y2 dx dy Q{y - x)3v a+1 p a (x)p a+2 (y), (2.5) 

where, as in the microscopic case, we have used the constraints (|2.2|) to rewrite 
the energy in a form analogous to (|1.6p . The constant term in (\2A\i , equal to 
— £ a r Q (l + r a+2 — r Q+ i)/2, is independent of the profiles p a (x) and has been 
added in to simplify the expression ()2.5|) . The free energy (multiplied by (3) 
J~({p}) associated to the canonical ensemble measure (jl.4p . is then given by 

F{{p}) = P£({P}) ~ S({ P }) 

= spy2 dx d y ®(y ~ x)v a+1 p a (x)p a+2 (y) 

a J ° J ° ( 2 - 6 ) 

+ T2 dx p a (x)\ogp a (x). 

a JO 

The functional J-({p}) is, up to an additive constant, the large deviation 
functional giving probabilities in the N — > oo limit [U El [TTJ, that is, the prob- 
ability of the profile p a (x) is proportional to exp [— ./V {F{{p}) — mm p T({p}))]. 
The typical equilibrium density profiles for the macroscopic model are thus those 
that minimize the free energy. There will be a coexistence of phases when the 
minimizer is not unique. To study this we have to consider solutions of the 
Euler-Lagrange equations (ELE) associated to the stationary points of (|2.6p . 

Let J- a (x) = 5 J- 1 '5p a (x) be the variational derivative of J- taken as though 
all of the p a are independent. Applying the constraints (|2.2p leads to the ELE 



Fa — Fc = constant, Fb — Fc = constant, (2-7) 
where the F a from (I2.6P are 

F a (x) = 1 + logp a (x) + 3/3v a+ ir a+2 

f x (2 8) 

-3/3 / dy(v a+1 p a+2 (y) - v a+2 p a +i(y))- 



5 



It is easy to show that these J- a satisfy 

a -p 

Y^ Pa (x)^(x) = 0, (2.9) 

a 

which with (|2.7p implies that J- a {x) is constant for all a. Then from f|2.8|) the 
ELE are given by 

p a (x) = p a (0) exp ^3/3 J dy (v a+1 p a+2 (y) - v a+2 p a +i{y))\ , (2.10) 
which may be evaluated at x = 1 to yield the boundary conditions 

Equivalently, one may also write the ELE in differential form as 
dp a 



dx 



3(3p a (v a +lPa+2 - V a +2Pa+l) • ( 2 - 12 ) 



The ELE are to be solved subject to (|2.2[) . Solutions of (|2.12p are stationary 
solutions of the hydrodynamic equations associated to the microscopic evolution 
of the model under the dynamics discussed in Appendix [A] [U [TT]; these hydro- 
dynamic equations have the form 

^ ~ 3 f Pa^) ■ (2.13) 



dt dx \ dx 

Examining the ELE (I2.12P one finds that there exist two constants of the 
motion, 

J>a(aO = l, (2-14) 

a 

and 

K = l[p^(x). (2.15) 

a 

To study solutions of the ELE we eliminate pc using (12. 14ft , so that the solutions 
give trajectories in the pa-Pb plane, and more specifically in the triangle 

PA,PB>0, PA + PB<1- (2.16) 

Any such trajectory lies within a level set of K; as we will see below, when the 
v a lie in region I such a level set is either the single point p a = v a or a simple 
closed curve encircling this point. When the v a lie in regions II and III the 
level curve is an open curve joining two vertices of the triangle. A solution of 
the ELE satisfying the boundary conditions ()2.11[) . or equivalently the constraint 
Jq Pa(x) dx = r a , is obtained by choosing first a value of K and then a portion 
of the trajectory labeled by K which is traversed in time one. See Figures [2] and 
[3] for typical level sets and trajectories. 
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In the special case where r a = v a the macroscopic free energy is rotation 
invariant (see Appendix [X] for details on the analogous result for the microscopic 
system). To verify this we consider rotated profiles 

[ p a (x — z) if x > z 

x = { ,,,: u - f " ■ ( 2 - 17 ) 

I p a \X + (1 — Z)) it X < Z 

The entropy is clearly unchanged by the rotation, S({p}) = S({p}), while the 
difference in energy is 



-£({/>}) = 3V / dx dy v a+ i (p a (x)p a+2 (y) - p a (y) p a+2 {x)) 

= dxv a+ i(p a (x)r a+2 -r a p a +2(x)) =0, (2.18) 

where we have used r a = v a and f^dx p a (x) = r a . As we shall see this case, 
which generalizes the r a = 1/3 case in the standard ABC model, plays a special 
role in the phase diagram. 

It was proven in [TJ Section 10] for the standard ABC model that solutions of 
the ELE always exist, and that the minimizer of the free energy must be given by 
one of the solutions. The same result holds for the generalized ABC model, with 
no modification of the proof required. The fundamental question is then whether 
or not there is a unique minimizer. This question was answered completely for 
the standard ABC model in [I], and we believe that the direct generalization of 
the behavior in that case should hold for the general model. We state this here 
as a conjecture, although it is partially established in [1] and in the remainder of 
this paper; the statement depends upon a critical temperature T c = P^ 1 , where 

& = 5 ; ^ ■ ( 2 - 19 ) 

Conjecture 2.1. Solutions of the ELE. (a) If the v a lie in region I and 
f a — v a for a = A,B,C then there exist (i) the constant solution p a (x) = v a , 

< x < 1, (ii) for (3 > n (3 C) n = 1,2, . . ., a unique solution corresponding to a 
trajectory which traverses one of the level sets of K exactly n times, which we 
refer to as a type n solution, and (Hi) no other solutions. The minimizer of the 
free energy is, for f3 < (3 C , the (unique) constant solution and, for (3 > f3 c , the type 

1 solution. At f3 c there is a second order phase transition from the homogeneous 
phase to the phase segregated, heterogeneous phase. 

(b) For values of the v a and r a other than those discussed in (a) there exists for 
every f3 a unique solution minimizing the free energy. 

Statement (a.i) here is a trivial observation, and the existence portion of (a.ii) 
follows immediately from the existence of a minimizer and the fact that f3 > (3 C 
the uniform solution is unstable (Section 14. 2ft ; we give an independent argument 
in Section [H Beyond this, as we will discuss below, we can prove all or part of 
this conjecture for some special values of the parameters besides the standard 
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case, v a = 1/3 for all a, for which the conjecture has been proven in full. In 
particular we can prove uniqueness of the solution of the ELE when (3 < 47r/3. 
This follows from the result in [I] that for such /3 the standard ABC model free 
energy functional is globally convex on the space of density functions satisfying 
(|2.2p . As the addition of external fields only adds terms that are linear in the 
particle densities to the standard ABC free energy (see (|2.5p ). the second variation 
of (|2.6p with respect to the density functions is the same as that of the standard 
model. Thus for (3 < 4-7r/3 the free energy for the generalized ABC model with 
external fields is also globally convex, implying that there is a unique solution 
of the ELE, which must be the minimizer of the free energy. For r a = v a this 
is just the constant solution p a (x) = v a , < x < 1. For other values of the r a 
there is for f3 < 4tt/3 a unique segment of a unique K = K(f3\r) trajectory which 
minimizes J-. 

For (3 > 4-7r/3 and r a = v a (so that the v a are in region I), we have the 
following additional results: 

(i) Conjecture 12.1( a) is proven when, for some a, v a = 1/2 and v a+ i = v a+ 2 = 
1/4, and checked numerically in other cases.. 

(ii) The constant solution is linearly stable for (3 < (3 C and unstable for (3 > (3 C . 

(iii) Nonconstant minimizers are always of type 1. 

(iv) For small enough K, or equivalently, for large enough /3, the type 1 solution 



When r a ^ v a for some a we prove uniqueness for the cases when one of the v a 
is zero and the other two have opposite signs, and when one of the v a is one. 

In the next section we describe general properties of the solutions of the ELE. 
The (with the v a in region I) is discussed in detail in Section HJ 

while the special cases when r a ^ v a are considered in Section [5j 

3 General properties of solutions of the ELE 

The trajectories of the densities p a {x) that are solutions of the ELE may be 
obtained by studying the level sets of the constant of the motion K (|2.15p . as 
described in Section [2l To do this let us define a line in the pa-Pb plane passing 
through the point (va, vb) by setting 



with m an arbitrary constant. The change in log A' as pa is varied along the line 
(13.11) can be manipulated using (|2.2p to yield 



is unique. 



PB = v B + m (p A - v A ) 



(3.1) 



d log if 
dpA 




(3.2) 



S 



Pa 



X 



Figure 2: Plots of pa (solid), ps (dashed), and pc (dotted), at right, and their 
corresponding trajectories in the pa-Pb plane, at left. See description under 
Region I, Case 1. 

Thus K is monotone increasing for pa < va and decreasing for pa > va- A 
similar result holds for ps and vb- The shape of the level sets of K, and thus 
also the trajectories of the densities, depends upon which region the v a lie in; 
compare Figures [2] and EJ We now give more details. 

Region I 

It follows from (|3,2p that in region I K achieves its maximum value at the 
center point (va, vb), where 

K\ pa=Va =K max = Hvl°, (3.3) 

a 

and decreases monotonically as one moves along any straight line in the (pa, Pb) 
plane starting from the center. K approaches its minimum value of zero on the 
boundaries of the triangle (|2.16p , where one or more of the particle densities goes 
to zero. As K is continuous in pa, Pb inside the triangle, this implies that the 
level sets of K in this case consist of a single point at the center (va, vb) and 
closed curves surrounding the center point. 

As the level curves of K are closed, nonconstant solutions p a (x) of the ELE 
must be portions of periodic functions of x. When the period r is greater (less) 
than 1, the trajectory of the densities makes less (more) than one full rotation 
around the center point. We will refer to solutions with n — 1 < 1/r < n, for n 
an integer, as type n solutions. For example, a solution that does not make one 
full rotation would be labeled as type 1, while a solution making exactly three 
rotations around the center would be labeled as type 3. This is consistent with 



the terminology used in Conjecture 12.1( a). The constant solution of the ELE is 
not assigned a type. 

Trajectories of the particle densities satisfying the ELE for a given set of v a 
depend upon the choice of the r a . For the v a in region I there are two cases to 
consider: (1) r a = v a for all a, (2) r a / v a for some a. 

Case 1: r a = v a . In this case the model is rotation invariant, and one has 
Pa(l) — y°o(0) for all a. Here both constant and nonconstant solutions of the 
ELE are possible. The constant solution is given by p a (x) = v a , < x < 1, as 
noted above. As /9 a (l) = Pa(0), see f)2.11[) . nonconstant solutions must have an 
integer number of periods in the interval x £ [0, 1], corresponding to the number 
of times the trajectory orbits the center. Note that, as one moves along the 
interval in x, the maxima (and minima) of the particle densities proceed in cyclic 
order, that is, after species a reaches its maximum (minimum) density, the next 
species to achieve its maximum (minimum) density is a + 1. 

Example numerical solutions of the ELE in this case and their corresponding 
trajectories with va = 1/2, vb = 1/3, vq = 1/6, and r a = v a for all a, are shown 
in Figure [2j In this plot the inverse temperatures are /3i = 13, 02 = 15, and 
/?3 = 20, all larger than /3 C , which is 4ir for this choice of the v a . The trajectories 
lie along lines of constant K, with K\ ~ .349, K2 ~ .291, and K3 ~ .189 for 
the solutions at /3i, /?2, and $3 respectively. Arrows indicate the flow along the 
trajectory as x increases, and the point (va, vb) is marked by a dot. For this case 
K max .364. 

Case 2: r a ^ v a . Here there is no rotation invariance and the densities at 
opposite ends of the interval are not the same, p Q (l) 7^ p a (0) for some a, see 
(|2.1ip . In this case only nonconstant solutions of the ELE are possible at finite 
temperatures, i.e. for (3 > 0. These solutions will be portions of the periodic 
solutions described in Case 1. In contrast to the however, solutions 

of type n do not exist for arbitrarily large values of n; there is some cutoff 
"-max > 2, which depends on the r a and the v a , above which type n solutions, 
n > "max, do not exist. This is because the average value of each density p a 
around one full orbit of the center is v a , so as n becomes large the average 
density for the full profile r a is steadily driven towards v a . See Appendix ICl for 
more details. 

Regions II and III 

When the v a lie in regions II and III, the point (va, vb) lies outside the 
triangle (I2.16p . Thus by (13. 2^ the level sets of K inside the triangle cannot be 
closed curves, as when the v a are in region I. Let us consider the value of K along 
one of the boundaries of the triangle, where p a = for some a. If v a is less than 
zero, K will be infinite along this boundary. If v a is zero K will be finite, and if 
v a is larger than zero K will be zero. On the vertices where two boundaries meet, 
with K equal to zero along one and infinite along the other, the value of K at the 
vertex is not well defined, depending on the way in which the limit is taken. Thus 
when the v a lie in regions II and III, lines of constant K for K finite and nonzero 
will be curves terminating on the vertices of the triangle where boundaries along 
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Figure 3: Plots of pa (solid), ps (dashed), and pc (dotted), at right, and their 
corresponding trajectories in the pa-Pb plane, at left. See description under 
Regions II and III. 

which K is infinite and zero meet. When the v a lie on the boundaries between 
different regions, e.g. for the special values of the v a considered in Section El the 
level sets of K will be curves with ends terminating on either the edges or the 
vertices of the triangle. Solutions of the ELE are therefore nonconstant at finite 
temperatures and are not periodic. 

Additionally, one may easily see by considering the ELE (|2.12p that when one 
or two of the v a are negative, the density of one particle species will monotonically 
increase in x, and another species will monotonically decrease. If for example 
va < while vb,vc > 0, ps will be monotonically increasing and pc will be 
monotonically decreasing. Note that if the v a lie in region II then as x increases 
from to 1 the maxima and minima of the particle densities proceed in cyclic 
order, as for region I systems, but when the v a lie in region III the order is 
reversed. That is, in region III a maximum (minimum) density of species a is 
followed by a maximum (minimum) of species a + 2. 

Example numerical solutions and the corresponding trajectory of the densities 
with va = 4/5, vb = 2/3, vc = —7/15, = rs = rc = 1/3, are plotted 
in Figure El In this figure f$i = 1, 02 = 3, and /?3 = 5. The trajectories lie 
along lines of constant K, with K\ ~ .306, K2 ~ .167, and K3 ~ .066 for the 
solutions at f3± , $2, and respectively. Portions of the level curves not traversed 
by the solutions are marked with dashed lines. Arrows indicate the flow along 
the trajectory, and the point (va,vb) is marked by a dot. 
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4 The case r a = v a 

If r a = v a for all a (so that we are necessarily in region I) then as already noted 
the constant solution p a (x) = v a , < x < 1, is always a solution of the ELE, 
and other solutions are given by traversing, exactly n times for a type n solution, 
one of the simple closed curves which is a level set of K. We show in Section \4. 21 
that the constant solution cannot minimize the free energy when /3 > /3 C , and 
in Section 14.31 that no type n solution with n > 2 can minimize the free energy. 
Thus the existence of a unique minimizer as described in Conj ecture 12 . 1 ( a) would 
be established if we could show that no non-constant solution can exist for (3 < (3 C 
and that for f3 > f3 c there is a unique type 1 solution, exactly what is proven in 
PQ for v a = 1/3 for all a. 

Now under a change of scale t = ftx the ELE fl2. 12f) become 

= 3Pa (v a+ ip a+2 - V a+2 p a +l) , (4.1) 

and type 1 solutions of (|2.12p correspond to solutions of (|4.ip which have pe- 
riod /3, the inverse of the temperature. The period t(K) of the solution of (|4.ip 
is easily seen to be a continuous function of K for < K < K max ; a per- 
turbative calculation as in pQ shows that lim^^<^ max t(K) = (3 C as K and we 
establish in Section [4.41 that \\mx\o T {K) = oo. The existence portion of Con- 
jecture I2.1( a.ii) follows immediately from these observations. If one could show 
that the period of the solution of (14. ip is a monotonically decreasing function of 
K for < K < i^max then uniqueness would be established. However, one may 
establish a correspondence between the family of systems (14. ip parameterized by 
the v a in region I and a large subset of the generalized Lotka-Volterra family of 
quadratic centers, specifically, those parameterized by (b, c) with b > in the 
parameterization of |13j (see Appendix [B] for details), and monotonicity of the 
period as a function of orbit size is for this family an open question [13] . In terms 
of the parameters used in (|4.ip . monotonicity was established in pQ when all v a 
are 1/3, and below we establish monotonicity when, for some a, v a = 1/2 and 
= f Q +2 = 1/4. Finally, monotonicity at sufficiently low temperature (or 
equivalently, small K) follows, for any v a , from the correspondence established 
in Appendix IBI and the results of [T3] . 

4.1 v a = 1/2, v a+1 = v a+2 = 1/4 

When the v a lie in region I and two of them are equal, say va = vb = 7 with 
< 7 < 1/2, one may find explicitly an autonomous evolution equation for the 
third. We begin with the expression 

K = PA {xV PB {xV Pc {x) l - 2 \ (4.2) 
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which together with (|2.2p may be used to obtain the densities pa and in terms 
of pc- 

p A (x) = \U~ Pc{x) ± ^(l-pc(x)) 2 - AKV-rpdx)-^ , (4.3) 

PB (x) = \ (l - Pc (x) =F ^(l-pc(x)) 2 -4K 1 /7 pc ,( x )-A^ . (4.4) 

where A = (1 — 2 7 ) / 7 . With (|4.3j) the reparameterized ELE (|4.ip for is 

Pc(t) = 37Po(t)(pfl(t)-pA(t)) 

= 3 7 pc(t) (±^(1 " Pc(t)) 2 ~ 4K 1 /7p c ( x )-^ . 



(4.5) 



Squaring (j4.5j) we have 



^fc(t) 2 + tfe(M*)) = 0, (4.6) 

with 

This is the equation for a zero energy particle of mass 1 confined in a potential 
well. The zeros of the potential correspond to the turning points for the particle. 
Numerical calculations of the period, performed for many different values of 7, 
find that t(K) is a monotonically decreasing function of K. 

With 7 = 1/4 and thus A = 2 the potential becomes quartic, and an analytic 
calculation is possible. Here the zeros of Uk are 



\±P0, po = ] J\-2K\ (4.8) 

Note that in this case 

Kmax = l/2-\/2, so < po < 1/2. The period t{K) of the 
type 1 solution p(t) may now be directly calculated, yielding 

rl/2+p 9 

t(K) = / dp 



l/2-po yJ-2U K {p) 
( 



16/3 _ p 



Vl/2 — Po — Po 



7T 4poVl/2-^ 
2' 



y/l/2-pl-po)' 



(4.9) 



Here F(ir/2, •) is the complete elliptic integral the first kind. The period (|4.9p is a 
monotonically decreasing function of if, with (as expected) lim^ /-K max (t(K)) = 
8\/2vr/4 = /3 C and limx\o t(K) = oo. Thus for all /3 > /3 C there is a unique value 
of K for which t(K) = f3 and hence a unique type 1 solution of the ELE, and 
there is no type 1 solution for (3 < f3 c . 
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4.2 Linear stability of the constant solution for j3 < f3 c 



We will now consider the linear stability of the constant solution p a (x) = v a , 
< x < 1, via a computation similar to that of [TJ. Let us consider two bounded 
continuous functions (ft a ( x ), (ftB(x), satisfying 



/ dx 4>a{x) = dx </>s(x) = 0, 
Jo Jo 



dx 4>a(x) 

and perturb the constant solution as 

(p A ,PB,Pc) ->■ (VA + e<j>A,V B + €(f> B ,Vc - e ((f) A + 4>b)) 



(4.10) 



(4.11) 



for some fixed e very small. Under this perturbation all the terms in the free 
energy linear in e cancel. The order e 2 contribution to the entropy is 

(4.12) 



dx 



—4> 2 A (x) + —(f) B (x) + —(4> A (x) + 4>b(x)) 2 
va v b vc 



The energy due to interaction with the external fields is linear in the densities, 
so the e 2 term for the energy is due entirely to the asymmetric mean field ABC 
interaction, 



/3 / dx dy Q(y — x) 
Jo Jo 



4>A(x)((t>A(y) + 4)b(v)) + <t>B(x)<f>A(y) 
-(4>a{x) + 4> B (x))4> B (y) 

= 3/3/ dx dy @(y - x)(ft B (x)(f) A (y)- 
Jo Jo 

Now expanding (ft a (x) and <ft B (x) in a Fourier series as 

oo 

(ft A (x) = ( a n sm (27rnx) + b n cos (27mx)), 

n=l 

oo 

4>b(x) = ( c ™ s ™ (2vnx) + d n cos (lirnx)), 



(4.13) 



(4.14) 
(4.15) 



»=i 



we see that the second variation of the free energy functional around the constant 
solution is given by 



£ 

n=l 



VA + V C , 2 . ,2x , V B +V C , 2 , ,2- 

(a n + b n ) + — — — (c n + d n/ 



vavc v B v c 
2 

H (OnCn + b n d n ) 

vc 



3/3 



7rn 



(4.16) 



If we write u n = (a n , b n , c n , d n ), the second variation can be expressed in matrix 
form as Y^Li u n Mu^/ (vav b v c ), where 



M 



v B (l - v B ) v A v B 3/3/(27rn) 

v B (l - v B ) -3/3/ (27m) ua^b 

ua^s -3/3/(27rn) - v A ) 

3/3/(27rn) ua^b - 



(4.17) 
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The eigenvalues of M are 



± 2 



va(1 - va) + v B (l - v B ) 



1 



(v A (l ~ va) + v B (l - v B )) 2 - 4 {vavbvc - 



(4.18) 



each with degeneracy 2. Thus for /3 < (5 C = 2tt/ (3^/vavbvc) the matrix M is 
positive definite, so the second variation of the free energy around the constant 
solution is also positive and the constant solution is a local minimum. At lower 
temperatures, j3 > /3 C , the smallest eigenvalue A_ becomes negative, so in this 
regime the constant solution can no longer be the minimizer of the free energy. 



4.3 Solutions of type n > 2 

In this section we will show that solutions of the ELE for r a = v a with two or 
more full periods in x G [0, 1] can never minimize the free energy. 

Consider a set of profiles {p}, not necessarily solutions of the ELE but satisfy- 
ing the constraints (|2.2|) . Then for any integer n > 2 we define the set of profiles 
{Pn} by 

j j + 1 

Pn,a{x) = p a (nx - j) for - < X < , j = 0, . . . , n - 1. (4.19) 

n n 

That is, we obtain the p n>a by shrinking the p a horizontally by a factor of n and 
repeating these reduced profiles n times in the interval [0, 1]. We claim then that 

S({Pn}) = S({p}), (4.20) 

£({p n }) = -8 ({p}) + (l - -) £ (W), (4.21) 
n \ n J 

where £({r}) = 9vavbvc/2 is the energy of the constant solution p a (x) = v a . 

The proof proceeds by direct calculation. First let us consider the entropy of 
the new profiles, 



S({Pn}) = ~2^2^ / p nta (x) log p n>a {x) 

a j=0 J i/ n 
n-1 j „! 

= "EE-/ dx Pa (x)log Pa (x)=S({p}). (4.22) 
Here we have used a change of variables with x = (x + j)/n. The calculation of 
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the energy proceeds similarly. Using the same change of variables we have 




dy Q(y - x) v a+1 p a {x)p a+2 (y) 



dy v a+ ip a (x)p a+2 {y) 



y+k-x-j 



n 



) 



v a+ ip a (x)p a+2 (y) 



(4.23) 



as there are n(n — 1) /2 pairs of indices with 0<j<k<n — 1. 

Now there are two cases to consider. If £({p}) < £({r}), then by (14.231) 
£({p n }) > £({p}), thus the original profiles p a have a lower free energy and 
thus the p n ^ a cannot be minimizers. If instead £({p}) > £({r}), then £({p n }) > 
£({r}) as well, and as long as p a is not the constant solution S({r}) > S({p n }). 
Then ^({r}) < F{{p n }), so again the p n ^ a are not minimizers. A type n solution 
of the ELE at an inverse temperature f3 is of the form p njCO where p a is a type 1 
solution at an inverse temperature f3/n. Thus no type n solution for n > 2 can 
minimize the free energy. 

4.4 K-f3 relation 

It follows from the result of Section 14.31 that when r a = v a a minimizer of the 
free energy must be either the constant solution of the ELE or a type 1 solution. 
Let K(f3) denote the value of K for the minimizer at temperature (5\ if several 
minimizers exist then we choose one of them arbitrarily to define K{(3). Of 
course, if K(f3) < K max then t(K{J3)) = f3. If t(K) is monotonic then K((3) 
will be the inverse of the function t{K) and must be continuous, since t(K) 
is, but we cannot show this and thus cannot rule out the possibility that K(f3) 
may be discontinuous. However, we do show here that K((3) must be monotonic 
decreasing. 

We begin with the integral form (|2.10j) of the ELE: 



log p a {x) = logp a (0) + 3/3 dy (v a+1 p a+2 {y) - v a+2 p a+1 (y)). (4.24) 
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Substituting this into the entropy (|2.3p and using p a (x) dx = r a gives 
Y^r a log/3 a (0) = -3/3 S2 / dx dy p a (x) (v a+ ip a+2 (y) - v a+2 p a+ i(y)) 

-s({ P }) 

= 6/3V]/ dx dy ®(y - x)v a+l p a {x)p a+2 (y) 
~ Jo Jo 

-S({p})-3PJ2 

VaTa+lTa+2 

a 

= 2{3£({p}) - S({p}) - 2(3£({r}), (4.25) 

where £({r}) (see (|2.6j) ) is the energy of the constant profile p a {x) = r a (which 
is not a solution of the ELE unless v a = v a ). (|4.25p is a general relation which 
holds for all profiles {p} satisfying the ELE, i.e. for all stationary points of J-"({p}), 
whether or not r a = v a . 

Now suppose that r a = v a for all a. The left hand side of (|4.25p is then just 
log-ff = ^2 a v a log p a (x), the same for all x £ [0, 1]. Then if {p} = {p^'} is the 
minimizing solution of the ELE corresponding to K(f3), (|4.25p becomes 

log #03) = 2P8({pW}) - S({pW}) - 2(3£({r}), (4.26) 

where now £({r}) = £({p^}) = 9vavbvc/2- For (3 2 > Pi we subtract the 
corresponding equations ()4.26p and rearrange terms to obtain 

log = 2(J^({pW»}) -TUiP^})) + {S({P^}) ~ S({ P ^})) 

+ 2(/3 2 - fk){8{{ p W>}) - £({p {0) })). (4.27) 

where we have indicated the explicit /3 dependence in (|2.6p by writing T$. But 
it follows from simple general thermodynamic arguments that both S({p^}) 
and £({p(^}) are monotonic decreasing functions of /3, and since p^ 2 ^ minimizes 
J-/3 2 all three terms on the right side of (|4.27p are nonpositive. This establishes 
the monotonicity of K((3). We note in addition that, if at a particular value 
of the temperature ft there exist several minimizers of the free energy, {pf^}, 
i = l,...,n, then log (Ki/Kj) = S{{pf ] }) - S{{pf ] }). 

One may also see that when f3 is large, K(/3) must be small. For from (|4.26l) 
we have 

= 2 (*({„(«}) - £({r})) + ^SUfW}); (4.28) 

since the energy and entropy are bounded functions and £({p^}) is decreasing 
in j3 the right hand side of (|4.28p approaches a finite value as (3 oo. In fact, one 
can show that lim^ £({p^}) is the ground state energy per particle 3vavbvc 
(see Appendix |A|) . so that asymptotically log K((3) ~ — (3/2)vavbvc(3- 
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5 Special cases when r a ^ v a , with the v a outside of 
region I 

For certain values of the r a and v a we are able to prove Conjecture 12.1( b). There 
are two such special cases. In the first case, one of the v a is zero and v a ±i have 
opposite signs (so that one of v a ±\ must be greater than one). In this case, the 
v a lie in fact on the boundary between regions II and III. In the second case one 
of the v a is one. We will present the proof of Conjecture 12.1( b) for these cases in 
the sections below. 

5.1 v a = 0, v a ±i > 1 

For definiteness let us take vc = 0, v B > 1. On the orbit with conserved quantity 
K we have p A (x) VA pb(x) 1 ~ Va = K and so 

p A (x) = K l '^p B {x)\ p B (x) = K- l l^-Vp A {x) l '\ (5.1) 

where 7 = {v A — 1)/^A satisfies 7 > 1. Thus the profiles satisfy 

p' A = 3f3p A (v B p c ) = 3/3(1 - v A ) PA (l - PA - K^-^p 1 / 1 ) = f K (p A ), (5.2) 
p' B = 3pp B (-v A p c ) = -3Pv APB (l -pB- K l ' VA p B ) = 9k {pb). (5.3) 

Note that if K* > K then 

fx(p) > fK*{p) > and g K *(p) > g K {p) > 0. (5.4) 

Now consider two profiles p a (x) and p%(x) satisfying the ELE which have 
different starting values: (p A (0), ps(0)) 7^ (p A (0), p* B (0)); we claim that the cor- 
responding averages are not equal: (r A ,r B ) 7^ (r* A ,r* B ). We denote the corre- 
sponding conserved quantities by K and K* and consider several cases. 

Case 1: K = K* . Without loss of generality we take p* A (0) > p A (0) and hence 
by (|5.ip also p* B (0) > Pb(0). Since p a (x) and p* a {x), a = A,B, satisfy the same 
differential equation which by (|5.4|) has a strictly negative right hand side we 
have p a {x) > p* a {x) for all x and so 

r A > r A and r* B > r B . (5.5) 

Case 2: p* A (0) = Pa(®)- Now without loss of generality we may take K* > K 
and so by (15. ip we have p* B (0) > Pb(0). Thus by (|5.3p and ()5.4p we have that 
p* B {x) > p B (x), and by (|5.2p and (|5.4p that p A (x) > p* A (x) for < x < 1, and so 

r* A < r A and r* B > r B . (5.6) 

The general case: Since we have dealt with the possibility K* = K in Case 1 
we may take K* > K. The case p^(0) = p A (0) has been considered in Case 2. 
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If p* A (0) > p A {0) we introduce a profile p** with p* A *(0) = p* A (0) and if** = K; 
then by Cases 1 and 2 we have 

r* B > r* B * > r B . (5.7) 

If instead p* A (0) < Pa(0) we argue similarly, introducing a profile p** with 
Pa (0) = Pa{0) and #** = K*, obtaining 

r A > r* A * > r* A . (5.8) 

Thus the solution of the ELE must be unique. 

The vq = 0, vb < 0, va = 1 — vb > 1 case may be argued very similarly, 
but making use of comparisons between the values of the densities at the end 
of the interval p a (l) rather than the initial values /5 Q (0), a = A, B. Otherwise 
the argument proceeds identically to the case above, so we will omit the full 
derivation here. 



5.2 v a = 1 

When one of the v a = 1 it is possible to solve the ELE exactly. Let us assume 
va = 1 and vb, vc are not zero, a trivial case, so the v a lie in region II. The ELE 
for pa and pb become 



dpA 
dx 

dpB 
dx 



3f3v B Pa(pb + Pc) 

3(3v B p A (l-p A ), (5.9) 
3Ppb(v c pa - Pc) 

-3P PB (1- Pb) + 3/3(1 -v b )paPb- (5.10) 



The equation for pa may easily be solved to obtain 

= l + C J-3/h*" ( 5 - U ) 

where ca = 1/pa(0) — 1. Using (|5.11j) one may then solve for the density ps, 



1 - v B 



1 



Pb(x) = c a (c A + e 3 ^*) VB - TJ — — , (5.12) 

with the constant 

CB = (i + CA )V«*( — 1 M --) 
V(l + ca)pb{0) c A ) 

To show that the solution is unique, we must demonstrate that there is only 
one choice of the initial values p a (0) which will yield a particular set of average 
densities r a at a given (3. This may be seen directly from the density profiles 
(15. lip . (j5.12p . which depend simply on ca and cb- It is easy to see that 

< o, %M < o, (5.13) 
dcA dc B 
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and as dcA/dpA(0) < 0, dc B /dps(0) < 0, we have 



^ >0, 3^T>0- (5-14) 



d PA {0) dp B (0) 

Thus there is a unique solution of the ELE, given by (|5,lip . (|5.12p . for all f3. The 
average density of the profiles is given by 

I g c A _|_ g3^B 

r a = —t, — log , (5.15) 

3(3v B B l + e c A ' y ' 

r B = 1 ^log- -n , (5.16) 

r c = l-r A -r B . (5-17) 

In principle one may invert equations ()5. 15[) , ()5. 16[) to write the solutions in terms 
of the average densities ta and r B rather than pa(0) and p B (0)- 
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A Properties of the microscopic system 

Many of the properties of the canonical measure pp for the standard ABC model, 
described in [TJ Sections 1 and 2] , hold with relatively minor changes for the gener- 
alized model studied here. Below we discuss a few of the necessary modifications. 

As in the standard ABC model, there is a certain nearest neighbor exchange 
dynamics which satisfies detailed balance with respect to the canonical Gibbs 
measure pp. In this dynamics a particle of type a at site i and a particle of type 
7 at site i + 1 exchange places 07—7-70 with rate q a ^, 

( e -Wv a+ 2/N if 7 = Q + l 

II it 7 = a — 1 

where the v a are as in (II. 7p . When the £ a are not all zero (i.e. the v a are not all 
equal), the system no longer has cyclic symmetry in the particle types. 

In analogy with what happens in the standard ABC model when N a = N/3 
for all a, we have here a rotation invariant energy when v a N = N a for all a 
(this is clearly only possible when the v a lie in region I). That is, if one imagines 
connecting site to site 1 and then rotating the configuration rj, the rotation 
leaves E unchanged. A simple way to check the rotation invariance is to consider 
moving a particle of type a from the end of the interval at site to site 1, and 
translating all the other particles from sites i to i + The change in energy after 
this rotation is then 3 (t> Q +i A r Q+ 2 — v a+ 2 N a+ {) /N, which vanishes for v a N = 
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N a . Note that in this rotation invariant case the rates (| A. 1|) for a particle of type 
a at a site % and a particle of type 7 at % + 1 to exchange become 

f -WN a+2 /N* if 7 = a+ l 

9a7 = r (A.2) 

II if 7 = a — 1 

which satisfies the general condition on exchange rates derived in [3] for which 
detailed balance holds on the ring. 

We remark that the energy in this rotation invariant case may be constructed 
in a different way, beginning from the standard ABC model energy with £ a = 0. If 
one wishes to write down an ABC-like energy that is explicitly rotation invariant, 
one way to achieve this is to take the standard ABC energy and average over 
starting the sum at each point of the lattice, i.e. 

N ( N+k N+k \ 

%) = Jf E [jf E E E - 0*k(0*fe+2(j) , (A.3) 

k=l \ a i=l+k j=l+k J 

where one imagines the interval with periodic boundary conditions, such that site 
N + m refers to site m. This gives 

1 N N ( ■ - \ 

%) = n EE E ( ®0" - + V ^w^+ 2 ^')' ( A - 4 ) 

a i=l j=l ^ / 

where the interaction term depends upon the distance between sites i and j as 
well as their order on the line. One may easily rearrange this expression to find 



N N 



%) = jf E E E - 3 Va{i) r?a+2(j) ' (A - 5) 

a i=l j=l 



which is identical to (|1.6p up to a constant when v a N = N a . 

Now let us consider the ground states of the model. As in the standard ABC 
model, in the (3 — > 00 limit the particle species become phase separated, with the 
ground states consisting of macroscopic domains of pure A, B, and C particles. 
The arrangement of these domains may depend upon the v a as well as the number 
of each particle species A^. For values of v a in regions II and III, where one or 
two of the v a is negative, the ground state is completely determined by requiring 
that all nearest neighbor configurations be stable, i.e., that the energy may not be 
lowered by making a nearest neighbor exchange. Consider the possible orderings 
of nearest neighbor pairs of different particle species: AB, AC, BA, BC, CA, CB. 
When the v a do not all have the same sign, the energetically preferred alignment 
of two nearest neighbor particles of different species will be cyclic for some pairs 
and anti-cyclic for others. Thus it is not possible to have an arrangement of 
four or more domains where all nearest neighbor pairs are preferably aligned. 
Generally, when the v a lie in region II (III), if t> 7 < (> 0), the ground state 
is given by three domains arranged in cyclic (anti-cyclic) order, i.e. in a cyclic 
(anti-cyclic) permutation of ABC, with the domain of type 7 particles in the 
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middle. For example, if v a < while vb, vq > 0, the stable nearest neighbor 
configurations are AB, CB, and CA, thus the ground state must consist of three 
domains of particles cyclically arranged as CAB. 

When the v a lie in region I, the ground states of the ABC model with exter- 
nal fields are closely related to those of the standard ABC model. To show 
the correspondence, let us consider a rescaling of the energy (jl.6p . E{rj) = 
E(rj) I '(vaVb vc)- The rescaled energy of a sequence of domains, where the 
ith. domain consists of k a i particles of type a (and no particles of any other 
type), is up to a constant factor formally identical to the standard ABC energy 
of the same sequence of domains, but with k a ^/v a particles in each domain. 
Of course, k a ^/v a may not be an integer, which must be taken into account 
when considering the degeneracy of ground states. One may however apply the 
same analysis as in p] to determine the lowest energy configuration of domains. 
This ground state energy per site is given, up to a constant term (see (|1.6p ). by 
3min {va tb tc, vb re va, vq M t\b}- One then finds that the ground state of the 
generalized ABC model, with v a in region I and N a particles of type a, is ar- 
ranged in the same order and has the same symmetries as the ground state of the 
standard ABC model with N a /v a particles of type a. There are three different 
cases. 

1. If one of the N a /v a is larger than the other two then there exists a unique 
ground state with three domains arranged in cyclic order, with the parti- 
cles for which N a /v a is largest in the middle. For example if Na/va > 
Nb,c/vb,c, the ground state arrangement of domains is CAB, with the 
particles of type A in the middle domain. 

2. If two of the N a /v a are equal and larger than the term for the third species, 
the ground state is degenerate, consisting of three or four domains cycli- 
cally arranged with the particles for which N a /v a is smallest placed on the 
boundaries of the interval. If for instance Nb/vb = Nc/vc > Na/va, the 
ground state is Na + 1 degenerate, with the domains arranged as ABC, 
ABCA, or BCA; the type A particles may appear on either the left or the 
right sides of the interval. 

3. If all of the N a /v a are equal, the ground state is N degenerate, consisting 
of arbitrary rotations of the domains arranged in ABC cyclic order. In 
this case, as discussed above, the energy of any configuration 77 is invariant 
under rotation. 

B The Lotka-Volterra family of centers and ABC-like 
systems of ODEs 

In this section we will demonstrate that the ABC model with external fields is a 
member of the Lotka-Volterra family of ODE systems when < v a < 1 for all a, 
i.e. when the v a lie in region I. 
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First let us change variables from x to t = 3/3 yJv~A vb vc x, a slightly different 
rescaling than that used in (|4.ip which will be convenient in the work that follows. 
With this change the ELE become 



Pa — Pa (u a +l Pa+2 ~ Ua+2 Pa+l) 

where u a = v a / ' ^Jva vb vq\ the u a satisfy 

UA + UB + UC = UAUBUC- 



(B.l) 



(B.2) 



After eliminating pc via (|2.14|) and uq via uc = {ua + ub)/(uaub — 1), we 
obtain the equations 



PA = PA 



PB 



-PB 



U A (1 + u|) 

U B (1 - pa) zr- PB 

UAUB ~ 1 

u B {l + u 2 A ) 



UA (1 - ps) 



UAUB - 1 



PA 



(B.3) 



In this form ua and ub are arbitrary positive constants satisfying uaub > 1, 
and the stationary point is 



va 



UAUB - 1 

ub (ua + ub) 



vb 



UAU B - 1 

ua (ua + «b) 



(B.4) 



We next shift the origin in the phase plane to the stationary point (|B.4p and 
make a linear change of variables: 



X 


= T 


'pa 


- VA 


JJ_ 




PB 


~ V B 



1 



UAUB - I 

In these new variables the equations become 



u B (l + u A ) u A (u A u B -1) 
ua (ua + ub) 



-y — bx 2 — cxy + by 2 , 



y = x + xy, 



where 



uaub - I 
1 + u 2 ' 



2ua + ub - u A u B 
1 + < 



(B.5) 
(B.6) 

(B.7) 



This is the canonical form for the generalized Lotka-Volterra family of centers 
presented in equation (1) of [T3] , 

Given parameters b and c with b > the equations (|B.7|) can be solved to find 



ua 



-c + Vc 2 + 46(6 + 1) 
26 



UB 



46(6 + 1) 



-0+^ + 46(6 + 1) 



c. 



(B.8) 



Clearly > and one finds easily that uaub > 1, so «b > 0. Thus the set 
of systems (IB.ip with ua, ub, > corresponds exactly with the family of all 
generalized Lotka-Volterra systems (|B.6P with 6 > 0. 
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Figure 4: Restriction on the existence of solutions of type n = 2 in the r a ^ v a 
case for va = 1/2, vb = 1/3, vc = 1/6, using (jC.ip . For values of (r^, tb) 
outside the bounded region, only type 1 solutions are possible. For n > 2, the 
bounded region in which type n solutions are possible would be smaller. 

C Restriction on the type of solutions for r a ^ v a , 
with v a in region I 



A naive estimate of the cutoff n max , as described in Section [3j may be made in the 
following way. Let us begin with a set of type 1 profiles, with average densities 
f a . These profiles will be a portion of the type 1 solution for the r a = v a case, 
stretched such that less than one full period of the p a (x) fits inside the interval in 
x G [0, 1]. We will define the length of the profiles Z, < Z < 1, to be the fraction 
of one full period of the p a {x) inside the interval. From these type 1 profiles we 
may make a set of type n profiles that satisfy the same boundary conditions as 
the original by rescaling x for the p a (x), such that the original type 1 profiles 
plus 7i — l full periods now appear in the interval. Then the average value of the 
densities for the new type n profiles, r n m will be given by 



f a l + nv a 
l + n 



(C.l) 



As a type n solution of the ELE for r a ^ v a will be a profile of this form, if the Z 
and r a in (jC.ip cannot be chosen such that the r n ^ a are equal to the specified r a , 
then this implies that a type n solution does not exist for that case. In Figure H] 
we show an example of the restriction imposed by this simple estimate. 
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